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A number-conserving hybrid Bose-Fermi variational theory 
is developed and applied to investigation of the BEC-Tonks 
gas crossover in toroidal and long cylindrical traps of high 
aspect ratio, where strong many-body correlations and con- 
densate depletion occur. 

03.75.Fi,03.75.-b,05.30.Jp 



The many-body ground state of a trapped atomic va- 
por Bose-Einstein condensate (BEG) is described in first 
approximation by Gross-Pitaevskii (GP) theory p],p| and 
in the next approximation by Bogoliubov theory|3] or 
time-dependent GP theory 0| . These fail if the conden- 
sate is appreciably depleted, as is the case near the BEG 
condensation temperature, or even at T = for suffi- 
ciently thin wave guides, low densities, and large scatter- 
ing lengths An extreme limiting case of this, the 
"Tonks gas" where transverse excitations are frozen, the 
dynamics reduces to one-dimensional (ID) motion, BEG 
disappears, and the occupation Nq of the lowest orbital 
behaves like with p < 1, can be treated exactly by the 
Fermi-Bose mapping method 0,^. The behavior of the 
BEG- Tonks crossover is of considerable interest since ex- 
periments are now approaching the Tonks regime [p|-p^ . 
An approximate theory in a strictly ID model, the Lieb- 
Liniger (LL) delta-function Bose gas ||l^ in a harmonic 
trap, has been given recently by Dunjko, Lorent, and 01- 
shanii , but it is unclear how to extend their approach 
to 3D as it is rooted in the ID LL model. The condensed 
fraction Nq/N is expected to be small for narrow waveg- 
uides due to the above-described behavior of iVo. 
Variational approaches such as Hartree-Fock-Bogoliubov 
(HFB) theory or its forerunner, the variational pair the- 
ory or Girardeau-Arnowitt (GA) theory of many-boson 
systems suggest themselves here. Here we gener- 
alize the GA theory so as to treat the gradual onset of 
"fermionization" as the Tonks limit is approached and 
apply it to the theory of the BEG-Tonks crossover. We 
also illustrate how the theory can be extended to 3D by 
allowing for variable transverse confinement. 

Number- conserving pair theory: The most general 
Bose pairing state involves excitation of pairs to ar- 
bitrary excited (uncondensed) states (n, m) by re- 
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peated application of pair excitation operators a^jflQ 
to the completely condensed iV-particle state \N) — 
(A^!)^^/^(aJ)^|0). Here we employ a number-conserving 
formulation |flo-n8[ and define unitary condensate anni- 



hilation and creation operators /3o 



-1/2 

(A'o + 1) oq and 



1) which commute with each other and 
with all annihilation and creation operators a„ and ajj for 
noncondensed atoms (n ^ 0); here Nq = ajag. Then the 
number-conserving generalized pair creation operator is 
al^aly^P^. In the case of a sufficiently symmetrical trap 
the general pairing state can be reduced to a ±fc pairing 
form 

-^|iV) , F=lyA,alalJl (1) 



const, e 
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by a suitable choice of orbitals Uk where u_fe denotes the 
time reversal conjugate of Uk and is real and even. 
This state is the vacuum of number-conserving Bose and 
Fermi quasiparticle annihilation operators Pq|l7|: 



ik = (1 



Al) 



-^/^0lak + AkalJo) 



^ 0, (2) 



and the general pairing state can be written as |<i>o) = 
;7|iV) where Lf^^ikU = a^. 

Toroidal and long cylindrical geometries: A convenient 
geometry for discussing the crossover is a toroidal trap 
of high aspect ratio R = L/Iq where L is the toroid 
circumference and £o the transverse oscillator length 
ia = y^h/muJo with uiq the frequency of transverse os- 
cillations, assumed to be harmonic. The transverse trap 
potential is assumed to be symmetric about an axis con- 
sisting of a circle on which the trap potential is mini- 
mum. The longitudinal (circumferential) motion can be 
described by a ID coordinate x in terms of plane-wave 
orbitals 4>k{x) = x~^/^e*'^^ satisfying periodic bound- 
ary conditions with periodicity length L, with allowed 
longitudinal quantum numbers kj = l-KnjjL with Uj = 
0, ±1, ±2, • • •. The corresponding 3D orbitals are taken 
to be Uki'r) = 4'k{x)(j}tr{p) where p is a transverse radial 
coordinate measured from the central circle of the toroid; 
these are cylindrical coordinates with cylinder axis bent 
into a circle of circumference L. This geometry can 
equally well be interpreted as an infinitely long, straight 
cylindrical waveguide with periodic boundary conditions 
in the longitudinal direction. The creation operators in 
the variational trial state |<i>o) refer to the 3D orbitals 
Uk{Y). Use of a single transverse orbital is justified both 
in the Tonks limit, where transverse excitations are ab- 
sent, and in the GP limit, where BEG is almost complete 
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and there is a single orbital determined by the GP equa- 
tion. 

Interatomic interaction: We use the usual Fermi pseu- 
dopotential f(ry) = 47ra(?i.^/m)5(ry) and assume that 
the s-wave scattering length a is positive. This leads 
to a well-defined problem in ID, the LL model [ p^ . 
Our toroidal system is "almost ID" since the transverse 
dimensions are confined, and we find that the varia- 
tional problem with the Fermi pseudopotential does not 
encounter the difficulties (divergences and poorly-posed 
variational problem) [ p^j20| found in the 3D case. A gap- 
less theory with more complicated pseudopotential, e.g. 

is not warranted here since we are concerned with 
the ground state, collective phonon excitations of the GA 
theory are gapless [^,^, and the quasiparticle gap re- 
moves an unphysical low-momentum divergence in the 
ID Bogoliubov theory. 

Pair Hamiltonian: The second-quantized many-boson 
Hamiltonian with units ?i = m = 1 is 



H 



k 



<^ka\ak 



qkk 



with single particle energy = (fc^/2) -I- e^^, transverse 
mode energy 
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where the transverse orbital is normalized according to 
/o \4'tr(p)\'^'^T^ pdp — 1, and interaction matrix element 
g = ina \4>trip)\'^'2'!Tpdp. The only interaction terms 
with nonzero expectation value in the state |<I>o) are 
those expressible in terms of momentum occupation num- 
ber operators Nk = a|afe and pair operators a-kdk and 
SfeaLfci and of those the transverse terms and the inter- 
action terms with q — sum to etrN + {g/2L)N{N — 1) 
where N — Nk is the total particle number operator. 
Replacing this operator by its eigenvalue N , one obtains 
a pair Hamiltonian JTot 

Hp = etrN + {g/2L)N{N - 1) 
+ 5][(fcV2) + {g/L)No]Nk 



k=ta 
9_ 
2L 



^ Y^miN, - l)Y'\plfa^kdk + h.c.} 



k' 



kk'^0 



+ (1 - ^kk' - ^/=,-fe')4«-fea-fc'"fe']' 



(5) 



whose expectation value in the state |$o) is identical with 
that of the full Hamiltonian (3). 

Bogoliubov theory: In the case of a an untrapped 3D 
system the Bogoliubov theory is valid at low densities. 
In the opposite limit of a strictly ID system with delta- 
function interaction (LL model), the Bogoliubov theory 



reproduces the leading terms in the exact ground state 
energy in the limit of high densities However, the 

Bogoliubov theory is not fully consistent cither in the 
strictly ID case or in our geometry: Wc find a Bogoli- 
ubov quasiparticle energy Uk = \k\\/nQg + and mo- 
mentum distribution function behaving like \2k\^^ y^griQ 
at low momentum, where no is the mean value of Nq/L. 
This leads to a logarithmic divergence of the depletion 
integral (fractional occupation of orbitals with |fc| > 0). 
We therefore proceed directly to the variational theory. 

Variational pair ( GA ) theory: The variational ground 
state energy Eqp is most easily evaluated by recalling 
that |<i>o) is the vacuum of number-conserving quasipar- 
ticle annihilation operators (2), applying Wick's theorem 
after rewriting Hp in terms of quasiparticle annihilation 
and creation operators via the inverse transformation 
(ilak = (l-A^)~i/2(^fe-Afe^Lj.). We pass to the thermo- 
dynamic limit in the longitudinal direction, i.e. N — > cx), 
L — !■ cx), N/ L = n (linear number density). Then the op- 
erator square root in the third line of (5) may be replaced 
by Nq with negligible error and Nq eliminated from the 
second and third fines via the identity Nq = N—J2k^o 
valid for eigenstates of total particle number with eigen- 
value TV. Applying Wick's theorem, passing to the ther- 
modynamic limit via X^fe^^o ~^ i^/'^'^) not- 
ing that all integrands are even functions of fc, one finds 
eventually 



eop — etr 
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> p 2 
[{^ + nog + -l2)Nk 



]dk, 



(6) 



where eop = Eop/N is the energy per particle, no = nf is 
the condensate number density with condensate fraction 
/ = wJ2k^o^k = I - ^dk, the momentum 
distribution function is Nk = ^\/{^ ~ ^'k)^ ^"^^ 



h 
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1 



dk 



Nkdk. (7) 



I2 = ng{l — f) and need not be evaluated separately. 
These equations are in one-one correspondence with Eqs. 
(21) and (22) of GA 10 via the correspondences — > 
0(fc), nog Pof^(fc), h ^ ^i(^), and I2 hik), but 
here the integrals Ii and I2 are independent of k due to 
our use of the Fermi pseudopotential. Minimizing the 
expression (6) by setting its functional derivative with 
respect to A^ equal to zero, taking into account the de- 
pendence of no, Ii, and I2 on A^:, one finds 

[nog - /i)(l + Al) - 2[{ky2) + nog + h]Ak = 0, (8) 

whose solution for A^ is 

Afe = [nog - Ii)-^[{k^/2) + n^g + h - Wfc], (9) 

where ujk is the quasiparticle energy 
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Substitution of this expression for back into the defi- 
nitions of Ii and I2 leads to integrals which can be eval- 
uated in closed form 
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dk fk^ 





(12) 



where A — Ii/n^g, and K = K{y/1 — A) and = 
E{y/1 — A) are complete elliptic integrals |^^. Here we 
have made use of the fact that /i > 0, since otherwise 
ujk would become imaginary at small k, which further 
requires that nog > Ji. The ground state energy per 
particle is finally found to be 



where 
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FIG. 1. Log-Log plot of the scaled energy per particle ver- 
sus the scaled linear number density n/g: Hybrid Bose-Fermi 
variational theory eos/g^ (solid line), exact LL solution 
eLh/g^ (long dash line), eop/g^ from the GA theory (short 
dash line). 

The short dash line in Fig. 1 shows a log-log plot of 
the scaled energy per particle eop/g^ as a function of 
scaled linear density n/g calculated using the GA the- 
ory in Eqs. (6)-(14) for constant coupling g. Fixing g 
corresponds to the regime where the transverse orbital is 
frozen as the lowest mode of the harmonic trap, in which 



(10) case our model coincides with that of Lieb and Liniger 
Jl^ . The long dash line in Fig. 1, almost indistinguish- 
able from the solid line, is the scaled energy per particle 
^LL{l)/g'^ = {"^1 gY^{l)l'^ from the LL theory. There 
is excellent agreement between the GA theory and LL 
theory in the high density regime n/g > 1, where the LL 

(11) ground state energy is known to agree with Bogoliubov 
theory but for low densities n/g <1 the predicted 
energies diverge. 

Hybrid Bose-Fermi variational theory: The failure of 
the GA theory for n/g < 1 can be traced to "fermion- 
ization" in which bosonic atoms become impenetrable at 
low densities [D; this is the ID Tonks gas limit where 
Bose-Fermi mapping applies. A variational method 
capable of treating "partial fermionization" is suggested 
to treat the BEC-Tonks crossover. In the crossover re- 
gion the system is intermediate between a BEG and a 
fermionized Tonks gas, so we model it as an interpene- 
trating mixture of a BEG of wN bosons with pair theory 
energy functional eop and (1 — w)N fermions with en- 

(13) ergy functional eoF- Assuming additivity of the Bose and 
Fermi energies as in the theory of ideal mixtures psf , one 
obtains an approximate energy functional (total energy 
per particle) 

eosin) ~ weop{wn) + {I - w)eoF((l - w)n) . (15) 

(14) The Fermi energy functional is just that of the ideal 
Fermi gas, eopin) = 7r^n^/6, since a contact interac- 
tion of spinless fermions is equivalent to no interaction. 
Then for a fixed linear density n one obtains eosin) < 
weop{wn) -|- (1 — w)"^7r^n^/6 and this can be numerically 
minimized with respect to w, using eop from Eq. (13). 
The solid line in Fig. 1 shows the resultant scaled ground 
state energy eop/g^ versus scaled linear density n/g for 
constant g, and excellent agreement is seen with the LL 
theory (long dash line) over the full range of densities. 
In particular, for low densities n/g the energy per parti- 
cle correctly approaches that of a free Fermi gas ^,|l^ . 
One notes from Fig. 1 that the variational energy lies 
above the exact LL energy, providing an a posteriori jus- 
tification of the minimization with respect to variation of 
w. 

BEC-Tonks crossover: Insofar as calculation of the 
energy via (15) is concerned, the system behaves as if 
(1 — 'w)N interacting bosons have been replaced by free 
fermions. The fermionic contribution dominates when 
n/g <^ 1 where the bosons are impenetrable and the 
Bose-Fermi mapping theorem applies, whereas the 
GA theory is accurate in the opposite limit n/g > 1. In 
the crossover region n/g 1 we define an effective Bose 
condensed fraction by wf = No/N where / = No/Nb 
and Nb = vuN. Figure ^ shows this condensed fraction 
as a function of the the scaled density for constant g, 
and a crossover occurs at n* /g = 2.17 where wf = 0.5: 
This crossover condition n*/g = 2.17 is very close to 
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Olshanii's prediction N* = L/TT\aio\ for the maximum 
number of atoms to form a ID Tonks gas, where we iden- 
tify N*/L = n* and g — 2/Tr\aiD\,aiD being the ID 
scattering length Petrov et al. derived the same 
condition for the crossover from a Tonks gas to a ID 
quasicondensate 

For an ideal Tonks gas w which requires n/g —^ 0, 
but as a measure of the density needed to approach a 
Tonks gas if we require 90% or more of the atoms to be 
fermionized (w=0.1) we need to satisfy n/g < 0.026. By 
comparison, for a condensed fraction greater than 90% 
(wf > 0.9) we require n/g > 30. By comparison, using 
hydrostatic equations applied to a trapped LL model , 
Dunjko et al. find that the density profile of the trapped 
gas is close to the Thomas- Fermi solution for n/g — 14 
(ry = 9) for which wf = 0.85, and close to that for a 
Tonks gas for n/g « 0.1 {rj — 0.07) for which wf = 0.25, 
where rj = n\aiD\ in their notation. There is therefore 
consistency between our results even though we do not 
consider a gas with longitudinal trapping. 




FIG. 2. Condensed fraction wf versus the logarithm of the 
scaled linear number density logio{n/g). 

In summary, we have developed a hybrid Bose-Fermi 
variational theory that accurately describes the BEC- 
Tonks crossover in ID. A key virtue of our approach is 
that it may be extended to 3D: Here we have assumed 
tight confinement and taken g — Aira \4)tr{p)\'^'2'^ pdp 
as a constant whose value for the unperturbed transverse 
ground orbital is 5 = 2a/£o where €0 = @- More 

generally, minimization of the ground state energy by 
variation of (j)tr subject to the normalization constraint 
leads to the following generalized GP equation: 



P4>tr 



__1 

+ 47ran(2 - 



ujlp^cjitr 



^ d \ ^ 1 

pdpj 2 

~2Xf + f^\^)\^tr?^tr, 



(16) 



with Eqs. (6)-(14) will allow for the study of the crossover 
from ID to 3D and will be the subject of future work. 
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